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STATIQUE DES FLUIDES 



I- Equations 

1-1 Equation integrate 

Soit un domaine D de fluide, de surface S, et au repos dans le repere d’etude R (x,y,z). En statique des fluides, 
on choisit en general la direction z, (de vecteur unitaire z ) confondue avec la verticale du lieu. Dans ce cas, 
l’orientation de z sera toujours « vers le haut » ( ascendante ). 




Mdesignera un point dans D : M e D ; n(M ) est le vecteur unitaire normal a S au point M. 
M* designera un point de S : M* e S ; est le vecteur unitaire normal a S au point M*. 

Le fluide exterieur a D est appele complementaire de D, et note D . 

Ce domaine fluide est soumis a : 



Des forces a distance (comme la gravite) 



df g (M) = p(M)F(M)d(o 



d-1) 



de densite massique F(M ) cxerccc sur chaque particule elementaire (au sens de la MMC) de masse p(M)d(o. 
Ces forces sont aussi appelees «forces de volume », c-a-d reparties sur tout le volume de fluide contenu dans D. 

Des forces de pression exercees par le fluide exterieur a D. Ces «forces de surface » sont exercees par le com- 
plementaire D de D en chaque point M* de la surface S et s’ ecrivent, en chaque point M* : 



df(M*) = dA 



( 1 - 2 ) 



ou on rappelle que n(M *) est la normale exterieure a D en M* c-a-d dirigee vers le fluide de D . 



La condition d’equilibre de D s’exprime a l’aide des deux torseurs d’actions exterieures : 



T(F/D) + T(D/D) = 0 



( 1 - 3 ) 



ou T(F / D) est le torseur des forces a distance, et, T(D / D) est le torseur des forces de pression. 
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Les elements de reduction de ces deux torseurs (resultante et moment) sont definis par : 



T(F / D)= { 

o 

T(D / D) = < 



R(F / D) = f p(M)F(M)da> 

JMeD 

M(0,F / D) = \om A P (M)F(M)da 

JMsD 

R(D/D)= J- p{M*)n{M*)dA 



M(0, D / D)= J 



= OM a- p(M*)n(M*)dA 



L’equation d’equilibre entre les resultantes s’ecrit : 



J p(M) F(M)dco+ 


{-p(M*)n(M*)dA = 0 


MgD 


M*eS 



1-2 Equation locale 

Decomposons le domaine D en « particules » elementaires parallelepipediques (« briques ») dont chaque face est 
perpendiculaire a l’une des directions du repere R, et dont le centre est le point M. Le volume d’une « particule » 
s’ecrit : 

da> = dx dy dz (1-5) 



Les forces exteri ern es appliquees a une « particule » sont les forces a distance, p (M)F(M) dx dy dz , et les for- 
ces de pression sur chaque face du parallelepipede dont nous cherchons les expressions ci-dessous. 

Si Lon appelle p(M) la pression en M, on a pour la face orientee par x et situee a l’abscisse x + dx : 



1 dp(M) 

p(M) + — dx 

2 ox 



x dydz 



et pour la face opposee situee a l’abscisse v : 
x dydz . 



1 dp(M) 

p{M ) r dx 



2 dx 

La force de pression suivant la direction x est la somme des deux expressions precedentes, soit : 
dp(M) 



dx 



- dx dy dz x 



( 1 - 6 ) 



On obtient le meme type d’expressions pour les quatre autres faces en permutant les vecteurs unitaires et les 
coordonnees. 



La condition d’equilibre s’ecrit done, pour la particule consideree : 

dx dy dz = 0 






dx 



dy 



dx 



( 1 - 7 ) 



c/oj = dxdydz ne pouvant etre physiquement nul, on doit avoir en tout point M de D: 



Mecanique_des_fluides_2A_MP.doc 



4 



21/11/2007 





p (M)F(M) 



0 



dp(M) _ dp(M) _ dp(M) _ 

x + y + z 

dx dy dz 



On appelle gradient de pression en M, le vecteur defini par : 



— A dp(M) dp(M) „ dp(M) _ 

grad p(M) = x + y + - ^ z 

ox dy dz 



D’ou l’equation locale (c-a-d au point M) : 



p(M )F(M)~ grad p(M) = 0 



( 1 - 8 ) 



(1-9) 



Si p (M) est constante en tout point Mde_D, ou si o(M) ne depend que de la pression , l’equilibre du fluide n’est 
possible que si F(M ) derive d'un potentiel scalaire U tel que : 



F(M) = -grad U(M) 



(I- 10) 



( U sera precise ulterieurement) d’ou l'equation fondamentale de la statique des fluides: 



-p(M)grad U(M) - grad p(M ) = grad [p(M) + p (M)U(M)]= 0 => p(M) + p(M)U(M) = Constante, VM e D 

( 1 - 11 ) 



Les surfaces equipotentielles ( U = Cte ) et isobares ( p = Cte) sont confondues. 

En chaque point M, la normale commune a la direction de F(M) . 

Remarque 1 : Resultante des forces de pression appliquees sur les faces d'un element de volume. Quelle 
que soit la situation (fluide au repos ou en mouvement) la somme vectorielle des forces de pression appliquees 
sur les 6 faces de l’element s’exprime avec le vecteur gradient conformement a l’equation 1-8: 



dF(M) = - 



grad p(M) 



dco 



Le gradient de pression est determine physiquement par l’equilibre des actions mecaniques existant dans le 
milieu fluide. Ainsi, en statique, il est egal au vecteur -p g dco qui n’est autre que le poids change de signe de 

l’element de fluide. En dynamique, le vecteur grad p(M) est couple a : (i) l’acceleration p dV / dt , (ii) la 

variation grad \p LT V/)] de l’energie potentielle du fluide en mouvement, (iii) aux forces de surfaces. 



Remarque 2: Equivalence entre l'integrale de surface et Fintegrale de volume. Sachant que 

(,- , 2) 



d’apres (1-9) (fluide incompressible), on a en reportant dans (1-4), 

J p(M) grad^^y dco+ ^ - p(M*)n(M*)dA =0 (1-13) 

MsD M*eS 



soit la relation remarquable : 



| p(M*) n(M*) dA = J grad p(M)da> 

M*sS MuD 



(1-14) 
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II- Fluide homogene incompressible dans un repere terrestre (galileen ap- 
proche). 

II-l Equation locale 

La densite massique de force a distance est alors F acceleration locale de la pesanteurg et derive du potentiel U 
tel que : 



g = -grad U = - 


~dU 5C/_ dU 

v + -^-y + z 

ox dy dz 


(II- 1) 


si z est la vertieale ascendante du lieu, les composantes du vecteur g dans R sont alors 


A 8U n DU dU 

0 = , 0 = , g = 

dx dy dz 


(II-2) 



On en deduit que U ne depend que de z et a pour expression : 



U(M) = gz(M) + U 0 



(H-3) 



oil U 0 est une constante d’integration. 

Puisque la masse volumique est constante en tout point de D, p(M) = p, V equation fondamentale SDF-3 
s’ecrit : 



p(M) + pgz(M) = Cte, \/M e D 



(11-4) 



On appelle pression motrice le terme p(M) + pgz(M) et on le note p „ (M ) . L’equation SDF-5 devient done : 



p g (M) = Cte, VM e D 



(H-5) 



Dans un fluide au repos, la pression motrice est constante. 



II-2 Equation des isobares et de la surface libre 

Sur une isobare, p(M) = C M =Cte, on peut done ecrire : pgz(M) = Cte , soit : z(M) = H M =cte, c’est 
l’equation d’un plan horizontal de cote H M . 

La surface libre est une isobare particuliere, sur laquelle p(M*)e st une constante par exemple la pression at- 
mospherique, et d’equation : z(M*) = H 0 . 
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II-3 Equation fondamentale integrate - Theoreme d’ARCHIMEDE 



Soit un corps K, de frontiere (surface) S et de volume 
V= J dco(M ) 

MeK 

immerge dans un liquide. 

D’apres le principe fondamental applique au domaine D de fluide qui remplace fictivement l’espace occupe par 
ce corps, on a : 

p g J da+ p(M*)d(M*)dA = 0 (II-7) 

M<eD 

D’autre part les forces de pression exercee par le fluide exterieur a D, (c’est-a-dire D ) sur la surface S 
s’ecrivent : 

R(D/K)= p(M*)Tj(M*)dA (II-8) 

M*sS 

On en deduit que : 

R(DIK) = -pjv\ (II-9) 

C’est la formule donnee par le theoreme d’Archimede : tout corps K immerge dans un liquide est soumis de la 
part de celui-ci a une poussee R(D / K ) verticale dirigee vers « le haut », de module egal au poids du volume V 
de liquide deplace par le corps. 

Interpretation. 

En statique, la presence d’un corps immerge ne modifie pas le champ de pression dans le liquide. Le champ de 
pression p(z) est identique, qu’il y ait ou qu’/7 n’y ait pas de corps. L’integrale II-8 est done identique en 
l’absence du corps, si on la calcule sur une surface (5) identique a cede que le corps occupera quand on va 
I’immerger. En l’absence du corps, la resultante R(D / K) est egale a l’oppose du poids du volume de liquide 
enferme dans (5) et en equilibre. D’ou 11-9 et la phrase consacree qui suit. 

III-Fluide homogene compressible dans un repere terrestre (galileen appro- 
che). 

III-l Equation locale 

On a toujours l’equation differentielle au point M : 
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(III-l) 



dp(M) 

= -p(M)g 



A J Cas cTun LIQUIDE compressible 

Soient : p(M) la masse volumique du liquide a la pression en un point M, 

Po=P( M o) et P( M o) 

la masse volumique et la pression en un point M a de reference dans le fluide suppose a temperature constante. 
L’equation d’etat du liquide s’ecrit alors en fonction du coefficient de compressibilite isotherme du fluide, % : 

P (M) = p (M 0 )[l + %{p(M) - p(M 0 ))] 

L’equation differentielle s’ecrit alors : 

P (M 0 )[l + x(p(M) - p(M 0 ))] g = - (HI-2) 

az 

elle s’integre par separation des variables : 

[l + y{p(M) - p(M 0 ))] g = e JpfM ° )g_ [ z(M ° } _ z(M) ] (III-3) 

ou z(M 0 ) et z(M ) sont respectivement les cotes de M 0 et de M. 



En general, l’exposant est tres petit devant 1, on peut alors remplacer l’exponentielle par son approximation au 
l cr ordre, d’oii : 



p(M) - p(M 0 ) = p (M 0 )g[z(M 0 ) - z(M)] 



t | Xp(M 0 )g z(M 0 )-z(M) 
2 



Si on choisit le point de reference sur la surface libre de cote z(M 0 ) = 0, et que M est a la profondeur 
h = -z(M) , on a : 



P (M) - p(M 0 ) = p(M 0 ) g h 



t | X P( M o )g h 
2 



(HI-4) 



Le terme de compressibilite ne devient significatif qu’aux profondeurs importantes. Pour l’eau, x = 5E -10 Pa 1 
dans les conditions normales, le terme XP^ 7 vallt 0,01 a partir de 2000 metres. 

B / Cas d’un gaz. 

Soient encore : p(M 0 ) , p(M 0 ) , T(M 0 ) , la masse volumique, la pression, la temperature, au point de reference 

M 0 ■ 

En considerant que le gaz est parfait, on a : 



p(M) = p(M 0 ) 



p(M) T(M 0 ) 
P(M 0 ) T(M) 



(HI-5) 



Dans un gaz a temperature constante, on a : 



Mecanique_des_fluides_2A_MP.doc 



21/11/2007 





P(M 0 ) 



pm 

P(M o) 



g = - 



dp(M) 

dz 




Cette formule permet d’estimer revolution (decroissante) de la pression avec l’altitude. 
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DYNAMIQUE DES FLUIDES 



I- Coordonnees de Lagrange et coordonnees d’Euler 

Les coordonnees de LAGRANGE sont un ensemble de trois parametres a x , a , , a 3 , qui caracterisent la posi- 
tion de la particule a l’instant initial t = 0. 11 s’agit des coordonnees dans le repere /ive de reference du point 
avec lequel coincide la particule a l’instant t = 0. 

Les coordonnees d’EULER sont un ensemble de trois parametres x, y, z, caracterisant la position de la particule 
a l’instant t considers 11 s’agit des coordonnees dans le meme repere fixe que precedemment, du point avec le- 
quel coincide la particule a l’instant t. 

On notera bien que la direction z du repere fixe n’est pas necessairement confondue avec la vertieale du 
lieu. 



Pour une particule donnee, done que Ton « suit dans son mouvement » , on a pour t > 0 : 



x = fi(a l ,a 2 ,a J ,t) 
y = f 2 (a i,a 2 ,a 3 ,t) 
z = f 3 (a u a 2 ,a 3 ,t) 



a, = g x (x,y,z,t) 
« 2 = g 2 (x,y,z,t) 
a 3 = g 3 (x,y,z,t) 



( 1 - 1 ) 



A l’instant initial, on a : 

fi(a l ,a 2 ,a 3 ,0) = a l 

f 2 (a 1 ,a 2 ,a 3 ,0) = a 2 (1-2) 

f 3 (a l ,a 2 ,a 3 ,0) = a 3 



II- Vitesse 



La vitesse d’une particule est decrite par ses trois composantes u, v, w, sur la base orthonormee (x, y, z) du re- 
pere de reference, e’est-a-dire : 




(II-l) 



(II-2) 



(H-3) 



(11-4) 
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Ill- Derivee « particulaire » 



Soit une grandeur physique scalaire (|> liee a la particule de fluide en mouvement (concentration, masse volumi- 
que, composante de vitesse,. . ..). 

Si cette grandeur est exprimee en fonction des coordonnees d’EULER de la particule (x,y,z), on peut ecrire : 



<KM,t) = <Kx,y,z,t) 



(HI-1) 



On appelle derivee particulaire de cj) , sa derivee totale par rapport au temps pour a, , a 2 , a } , donnes, soit : 



ck |> _ 3(j) 5(|) dx <5<|> dy S<|> dz 

dt dt dx dt dy dt dz dt 

= ^r + (V-gmd)^ 



(III-2) 



V est le vecteur vitesse de la particule situee enMa l’instant t. 

(V ■ grad) designe l’operateur de convection applique ici a la grandeur (|> . 

Dans le cas d’une grandeur vectorielle quelconque, A(M, t ) = A(x, y, z, t) , le vecteur dA/ dt a pour composantes 
les derivees particulaires de chaque composante : 



dA 

dt 




= -^ L + (V-gr-dd)A x 



dA y 
1 dt 



= -^-+(V -grad )A y 



dA_ 

dt 



n A — *■ 

^ + (F-gradK 



(HI-3) 



IV- Acceleration 



Lorsque la grandeur vectorielle A est la vitesse de la particule au point M a l’instant t, on obtient le vecteur ac- 
celeration et ses trois composantes cartesiennes : 





du du - : 




ww* (Vg,!,i) “ 


- dV dV - 


t 

> 


r= ^r = ^r +(F ' grad)F=< 


~dt = ~dF + " ' grad) v 




dw dw - * 




wir +(f sn ‘ d,w 



dV 

Le terme represente l’acceleration locale ou temporelle. 



Le terme (V ■ grad)F represente l’acceleration convective. 

Le vecteur acceleration peut se mettle sous la forme qui fait apparaitre la vorticite 




(IV-5) 
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(IV-6) 



dV 

dt 



dV 

dt 

dV_ 

dt 



+ grad 
+ grad 



V 2 — >_ 

— — I- (rot V) aV 

V 2 _ - 

— + 2waV 



V- Champ des vitesses dans un ecoulement permanent (ou stationnaire ). 

On dit que l’ecoulement est permanent ou stationnaire si ses proprietes ne dependent pas explicitement du temps. 

En particulier, la vitesse ne depend pas explicitement du temps, c’est-a-dire que les composantes sont des fonc- 
tions uniquement des coordonnees d’ EULER x, y, z. 

On a done : 



du _ dv _ dw 
dt dt dt 



(V-l) 



Les composantes du vecteur acceleration se reduisent a : 
r,=(K-grad)ii r y = (F-g^d)v T z = (V ■ grad) w (V-2) 

L’expression du vecteur acceleration ne contient alors que [’acceleration convective : 



dV — j 

nr grad - 



V 



+ (rot E) a F 



(V-3) 



Si 1’ ecoulement est stationnaire les trajectoires et les lignes de courant sont confondues. 



VI. Ecoulement irrotationnel. 



On dit qu’a un instant t donne, l’ecoulement est irrotationnel dans le domaine D si le rotationnel de la vitesse 
V (M, t ) est nul en chacun des points M de D. Le rotationnel est un vecteur (defini en cinematique des fluides) 
qui caracterise en un point M, la rotation dans l’espace de la particule fluide. 



11 s’ecrit symboliquement : 



rot V (M,t) = V aV(M ,t) 



(VI- 1) 



II s’ecrit par exemple en coordonnees cartesiennes : 



rot V(M,t) 



dw(M,t) 


dv(M,t) 




du(M,t ) 


dw(M,t ) 




dv(M,t ) 


du(M,t ) 


[ dy 


dz 


x + 


dz 


dx 


y + 


dx 


dy \ 



(VI-2) 



VII. Equation locale de conservation de la masse. 

Soit un domaine D de fluide que 1’ « on suit dans son mouvement ». 

La masse m de liquide contenue dans D est constante dans le temps et dans l’espace (c-a-d au corns de 
l’ecoulement). 
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La derivee particulaire de m est done nulle : = 0 

dt 

On ecrit m sous forme de Lintegale de volume : 



m = 



J p (M,t)dm 

MeD 



dm 

dt 



d_ 

dt 



J p(M,t)dai 

MeD 



5p( ^ ,0 +div[p (M,t) V(M,t)] 



d(i>= 0 



(VII- 1) 



Cette relation 1 s’applique a tout element de volume, et done en particulier au volume elementaire da> ; ainsi, en 
supprimant le signe d’integration, on obtient L equation dite equation de continuity ou de conservation de la 
masse : 



5p( ^’ 0 +div[p(M,f) V(M,t)]= 0 



(VII-2) 



Comme 



div 



= pdivF + 



(f • gradjp 



on peut faire apparaitre la derivee particulaire de la masse volumique, fonction scalaire du point M et du temps t : 



dp(M,t ) 
dt 



+ p(M,/)e = 0 



8 = div V(M,t) 



(VII-3) 



Lorsque la derivee particulaire de la masse volumique dpi dt est nulle en tout point, on a divF(M,/) = 0 et 
Lecoulement est isovolume. 

On montre en effet que e represente le taux de deformation cubique (ou vitesse de dilatation ) de la particule si- 
tuee au point M. 

II est important de se souvenir dans quel cas dpi dt = dp/ dt + (V ■ grad)p = 0 . 

Statique des Fluides : dp / dt = 0 et V = 0 => dp / dt = 0 ; 

Fluides homogenes, incompressibles et indilatables : dp/dt = 0 et gradp = 0 => dpi dt = 0 ; 

Lignes de courant confondues avec les lignes isochores (ou p = Cte ) et masse volumique constante dans le 
temps (lluide inhomogene) : dp/dt = 0 et V est perpendiculaire a grad p => dpi dt = 0 

Lorsque Lecoulement est permanent ou stationnaire (le mouvement ne depend pas explicitement du temps), la 
derivee partielle de la masse volumique par rapport au temps est nulle : 



^"■'Ldiv 

dt 


\p(M,t) F(M,t)] 


= div 


[p (M,t) V(M,t)] 


|=0 



(VII-4) 



i ^ ,, ,, - _ _ _ , du 8v dw 

On appelle « divergence » d un vecteur V = ux + vy + xvz la quantite (scalaire) : div V = — — + — + . 

dx dy dz 
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On dit alors que le champ de vecteurs p V est conservatif. 

VIII. Equation integrale de conservation de la masse. Debit de liquide a 
travers une section de passage. 

VIII-1 Equation integrale 

Reprenons 1’ equation de conservation : 

| 5p( ^’ ?) + d |v [p(M,0 L(M,f)] d(£>= ^ P( ^ ,t) dw+ Jdiv[p (M,t)V(M,t)}l<D (VIII-1) 

MsD'- M<eD MeD 

et transformons la deuxieme integrale de volume en une integrale de surface en exprimant que la divergence du 
vecteur p (M, t ) V ( M , t) dans le domaine D est egal au flux de ce vecteur a travers S: 

J div[p(M, t ) V(M, t)}/co = J p (M*, t) V(M*, t ) n(M*)dA (VIII-2) 

MsD M*eS 

Or la projection 

V(M*, t) n(M*)dA = V(M*, t)dA = dq r (VIII-3) 

est le debit elementaire volumique dq v a travers l’element de surface dA oriente par n(M*) . Ainsi on a : 

| p(M * ^ dq v (M*) = 0 

MbD M*eS 

Lorsque le mouvement du fluide (ecoulement) ne depend pas du temps, on a cp / at = 0 et done : 

(VIII-4) 

Dans les applications ou le fluide est incompressible, p = p = Cte peut etre « sortie » de l’integrale : 

(VIII-5) 

L’integrale de surface J dq v ( M *) est tout simplement le debit en volume total a travers l’enveloppe S du do- 

M»<sS 

maine fluide D. On voit que ce flux total est nul. 

Si la masse volumique n’est pas constante, il faut raisonner sur le debit en masse, dq m = p dq v , d’ou : 

(VIII-6) 
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VIII-2 Debit a travers la section droite d’un tuyau. 

Ici le domaine D de fluide est un tronijon de tube constitue par une surface amont SI, une surface aval S2, et une 
surface laterale SL. Pour le debit en volume, en tout point de SL la projection 
V(M*,t ) n(M*)dA = V(M*,t)dA = dq v est nulle, puisque le liquide ne peut traverser la paroi (la vitesse est per- 
pendiculaire a la normale). 11 en est de meme pour le debit en masse. 

Le debit en volume est utilise en pratique pour les fluides incompressibles. 

Le debit en masse est utilise en pratique pour les fluides compressibles (gaz). 




on a done : 



*7 V ~ VI + ^ V2 ~ ^ tfv\— tfv2 OH ’ 


dvi 


= 




= CTE 


Le debit en volume est constant a travers 




les sections droites 


C lm = dmX + dm2 =0 => <1 ml = m2 OU \ 


d,n\ 


= 


C lm2 


= CTE 


Le debit en masse est constant a travers 



les sections droites 



Remarque : le signe du debit depend de l’orientation choisie pour la normale a la section droite. Si n est dans le 
sens de l’ecoulement, le debit est positif, negatif dans le cas contraire. 

Si on raisonne (cas frequent) sur les valeurs absolues, on utilise la valeur absolue \q v \ pour le debit en volume, 
\q m I pour le debit en masse. 



IX. Principe fondamental de la dynamique applique a un fluide parfait- 
Equation d’EULER 

On cherche la relation entre le mouvement des particules 
On applique ci-dessous le PFD a un domaine D de fluide 
M est un point de D. 

S designera l’enveloppe (ou surface) de D. 

M* est un point de S 

Le mouvement est etudie dans un repere absolu R(0, x, y, z) 



fluides et les forces qui s’exercent sur ces particules. 
que Lon suit dans son mouvement. 
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IX-1 Theoreme de la resultante en fluide parfait 



A chaque instant, la quantite d’acceleration A /R (D, t ) du domaine D est egale a la somme des forces exterieures 
qui s’appliquent sur D (on reprend les notations vues en statique des fluides): 



A ir (D, t) = R(F i D) + R(D / D) 



(IX-1) 



Ces forces comportent les forces de gravite R(F / D) de densite massique g , 



R(F / D) - J p(M,t) g da 

MeD 



(IX-2) 



et la resultante R(D / D) des forces de surface exercees par le complementaire D de D (fluide exterieur) sur 
l’enveloppe S , 



R(D/D)= J- p(M*,t)n(M*)dA 

M*<=S 

ou est la normale exterieure (unitaire) au point M* de S. 



(IX-3) 



R(D / D ) represente la resultante des forces de pression qui s’exercent sur S a l’instant considere. 



La quantite d’acceleration d’une particule etant egale a p (M, ^ dco , la quantite d’acceleration de D (ou 

dt 

resultante des quantites d’acceleration du domaine D) est donnee par l’integrale : 



A, R p{M,t)= Jp (M,t) dVi ^’ t] do, 



(IX-4) 



Finalement, le PFD stipule que : 



J p (M,t) ^ _ j" p (M,t)g d<o+ J ' - p(M*,t) n(M*)dA 



(IX-4) 



L’integrale des pressions sur la surface peut s’ecrire aussi sous la forme d’une integrate de volume (ceci se tra- 
duit physiquement par la poussee d’Archimede en statique des fluides) : 



^ - p(M*,t) n(M*)dA= J-grad [p(M*,t) \dA 



(IX- 5) 



11 s’ensuit 



p (M,t) J) _p (M,t)g +grad p(M*,t ) n(M*)da> 

dt 



= 0 



Localement, pour une particule fluide de centre M, l’equation du mouvement s’ecrit done : 



P (M, t) dV = p (M, t)g- grad p(M, t) 
at 



(IX-6) 



Cette relation est appelee equation d'EULER. 

L ’acceleration de la pesanteur s’ecrit en introduisant le potentiel scalaire U(M): 
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(IX- 7) 



g = -grad[(7 (M)\ 
U(M) = g Z(M) 



ou Z (M) est la cote du point Msur l’axe vertical ascendant du lieu (et non la coordonnee z). 
Si la masse volumique est constante, on a : 



dV(M,t ) 
dt 



-grad[ U(M )] - grad 



p(M,t) 



= -grad 



U(M) + 



p(M,t) 
P (M,t) 



Utilisons l’expression de [’acceleration : 



dV 

dt 



dV — t 
= — + grad 



V 2 — - 

+ (rot V) aV 



(IX-8) 



dV(M,t) 

dt 



+ grad 



V\M,t ) 



+ (rot V(M,t)AV(M,t ) = -grad 



U(M) + 



p(M,t) 

p(M,t) 



8V{M,t) 

dt 



+ (rot V(M ,t) aV (M ,t) = -grad 



V 2 (M,t) 

2 



+ U(M) + 



p(M,t) 
P (M,t) 



Posons : 



H(M,t) 



V 2 (M ,t) 
2 



+ U (M ) + 



p(M,t) 

p(M,t) 



le premier membre de 1’ equation precedente verifie done : 



+ (rot V(M, t) a V(M, t ) = -grad H(M,t) 
dt 



(IX-9) 



(IX-10) 



IX-2 Integrale premiere de l’equation du mouvement stationnaire 
Equation de Daniel BERNOULLI (1700-1782).. 

En mouvement stationnaire, on a en tout point M : 



dV{M,t ) 
dt 



(IX- 11) 



De plus comme 



dH = grad H(M, t ) • dl 



(IX- 12) 



la fonction H doit verifier : 



- (rot V(M, t) a V(M, t)-dl = dH 



(IX- 13) 
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Ainsi, selon que la direction du vecteur deplacement elementaire dl est confondue avec la vitesse V(M,t ) oil 
le rotationnel rot V (M, t) , on obtient respectivement, pour les integrates H. les proprietes suivantes : 

di colineairea V(M,t) => H(M, t) = Constante sur une ligne de courant ; 

dl colineaire a rot V(M,t) => H(M, t) = Constante sur une ligne tourbillon 

Une ligne tourbillon est tangente en chacun de ses points au vecteur tourbillon. 

Dans un ecoulement rotationnel. He st differe d’une ligne tourbillon a une autre. 



IX-3 Equation de BERNOULLI 

On appelle equation de BERNOULLI l’integrale constante de l’equation du mouvement sur une ligne de cou- 
rant. Elle s’applique au cas d'un fluide parfait, en ecoulement permanent et irrotationnel. 

On a done pour une ligne de courant donnee Li: 



VM e Li, H(M, t) 



V 2 (M,t) 

2 



+ U(M) + 



p(M,t) 

p(M,t) 



= Constante 



(IX-14) 



Si, dans l’ecoulement, le rotationnel est nul en tout point , alors la fonction He st la meme pour toutes les lignes 
de courant, c’est-a-dire dans tout le fluide : 



X- Aspect energetique de l’equation de BERNOULLI 

Raisonnons sur le travail des forces subies par la particule de volume da> lorsqu’elle effectue un deplacement 
donne par le vecteur deplacement 



d l = dx x + dy y + dz z 



(X-l) 



Ces forces sont : 

-Les forces de pression au point M, d’intensite dF p 



(M) = 



-grad p(M) 




-Les forces de pesanteur : dF g ( M ) = -pi/co g z = pi/co g = -pgradU(M) dco . 



(X-2) 

(X-3) 



Lorsque la particule effectue un deplacement dl , le travail des forces de pression et de pesanteur est egal a la 
variation d’energie cinetique : 



8 W p +8W g =8 E c 



(X-4) 



8E C =d[pdcoV 2 (M)/2 



8 W p =dF p (M)-dl = 



-grad p(M ) 



dayd l 



8W = dF (M) -dt = p 



-grad U(M) 



davdll 



(X-5) 



On se rappellera ici que la differentielle d’une fonction scalaire /(,V/) peut s’ecrire : 



df(M) = grad /'( M) ■ dl 



(X-6) 



Le travail de la pression et de la pesanteur s’ecrivent done : 
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8W p = dp(M) dco 


8W g =pdU(M)da = d\pU(M)\d(o 


Soit : 




t/[pAo V 2 (M)1 2 ]= dp(M)da + d\pU(M)\](o 

0 = d[p(M) + pU(M) + p V 2 (M) / 2 ] do 


(X-7) 



Comme le volume da de la particule ne peut etre nul, on en deduit que c’est la differentielle qui est nulle : 



d[p{M) + pU(M) + p V 2 (M) / 2] = dPt = 0 



Pt(M) = p(M) + pU(M) + ^pV 2 (M) 



La differentielle d’une fonction constante etant nulle, on en deduit que la fonction Pt(M) , que Ton appelle 

pression totale, est constante : 



Pt(M) = p(M) + p U(M) + p V 2 (M) / 2 = Constante 



(X-8) 



C’est l’expression mathematique du theoreme de Daniel BERNOULLI : dans un ecoulement irrotationnel de 
fluide parfait incompressible, dans le champ de pesanteur, la pression totale est constante sur une ligne de cou- 
rant. 



Dans cette demonstration, on s’est attache a une particule fluide donnee. Le resultat concerne done implicitement 
la meme ligne de courant. 

Si l’ecoulement est de plus irrotationnel en tout point, la pression totale est identique en tout point. 



XI- L’equation de BERNOULLI, en pratique 

D’abord : 



On remarquera que les trois termes de /’ equation de Bernoulli doivent se compenser en tout 
point M de la ligne de courant consideree, de maniere a ce que la pression totale soit CONS- 



TANTE. 



Ensuite : 

—1 —2 

La pression totale a la dimension d’une pression ( ML T ), son unite SI est le Pa. 

La pression totale etant homogene a une pression, elle peut etre exprimee par une hauteur de fluide equivalente 
qui a la dimension d’une longueur (L). 

Selon l’unite de travail choisie, les differents termes de l’equation de Bernoulli sont denommes comme suit : 



—1 —2 

Formulation homogene a une pression ( ML T ) 


Formulation homogene a une longueur ('Ll ou 




hauteur equivalente 


Pt(M) = p(M) + pg Z(M) + pV 2 (M)l 2 


H(M)= pt ' M) +Z(M)+ V 




Pg 2 g 
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XII- Pression motrice dans un ecoulement unidirectionnel stationnaire 

On rappelle que la pression motrice est definie en un point M a l’instant t par : 
p g (M) = p(M) + pgZ(M)\ (XII-1) 

et que Z(M) est la cote du point M sur la verticale ascendante du lieu. 

11 n’y pas forcement une relation simple entre z(M) coordonnees de M sur l’axe z du repere d’etude et Z(M), 
contrairement a la statique des fluides. Dans la figure ci-dessous on a par exemple represente un ecoulement des- 
cendant. SZ / 8x representera la pente de l’ecoulement. 




Considerons le cas simple mais instructif d’un ecoulement stationnaire de fluide incompressible de la forme : 
V(M) = u(M)x (XII-2) 

Les vitesses sont toujours colineaires a x . L’equation locale de continuite se reduit ici a divF = 0 , d’ou : 

-* Qn 

div V = — = 0 => u est independante de x (XII-3) 

ox 

Ainsi dans un tel mouvement, la vitesse ne peut etre que constante ; c’est un ecoulement uniforme. 

Les composantes r\, et f r de l’acceleration sont nulles ; seules f A peut etre non nude : 
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(XII-4) 



dv 

dt 




-^| t +p * z ] =0 

^|[p + P«z]-0 



p + pg Z est independante de y 
p + pgZ est independante de z 



Ainsi la pression motrice p = p + pgZ est constante perpendiculairement au plan de l’ecoulement uniforme. On 

dit alors que la pression p a line repartition hydrostatique car cette repartition est identique a celle donnee par la 
statique des fluides. 
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XIII. Ecoulement d’un fluide visqueux newtonien incompressible. 

XIII-1 Ecoulement interne et ecoulement externe 



Les ecoulements de fluide visqueux peuvent etre divises en deux categories : 

1/ Les ecoulements internes qui sont limites par des parois (ecoulements dans les tuyaux, canaux, tuyeres. . .) ; 

2/ Les ecoulements externes, qui sont des mouvements de fluides autour de corps tels que des avions, fusees, 
voitures, etc... 



XIII-2 Ecoulement laminaire et ecoulements turbulent 

Un ecoulement turbulent peut etre caracterise par une orientation aleatoire (ou fluctuation) des vecteurs vitesses 
en chaque point ; autrement dit chaque composate u(M,t), v(M,t), w(M,t), de V(M,t) obeissent a des lois de dis- 
tributions (au sens probabiliste) en fonction de l’espace et du temps. On dit que le regime de l’ecoulement est 
turbulent. Un ecoulement turbulent peut posseder une coniposante moyenne U (le mouvement global se fait 
d’ouest en en est par exemple) pouvant dependre du temps mais de maniere plus « reguliere »; les fluctuations 
ont alors lieu autour de cette composante moyenne, de sorte que la vitesse s’ecrit : 

V (M, t) = U + v' (XIII-1) 

C’est alors sur v' que Lon fait apparaitre des lois statistiques. 

Dans un ecoulement laminaire, le mouvement du fluide s’effectue en « couches » parallele entre elles, chaque 
couche possedant sa propre vitesse ; le profil de vitesses dans le fluide en mouvement est bien ordonne. Les 
vecteurs vitesse conservent une orientation stable au cours du temps. On dit que le regime de l’ecoulement est 
laminaire. 

Lorsque le mouvement laminaire du fluide degenere en un ecoulement turbulent , il perd son caractere ordonne 
et stable ; on dit qu’il y a transition du regime laminaire vers le regime turbulent, ou plus simplement transition 
laminaire turbulent. 

La difference entre le regime turbulent et le regime laminaire peut etre representee par les courbes donnant la 
vitesse en fonction du temps en un point Mfixe. 



Evolution des vitesses avec le temps en un point fixe 




temps 



XIII-3 Ecoulement laminaire etabli entre deux plaques paralleles infinies 

Le fait de considerer un ecoulement stationnaire permet de se ramener a un probleme independant du temps. 
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Le fait de considerer un ecoulement etabli permet de se ramener a un probleme independant de l’abscisse x : le 
mouvement n’est ni accelere ni ralenti selon la position sur l’axe x ; u ne depend pas de x. 

Le fait de considerer des plaques infinies permet de se ramener a un probleme plan (a deux dimensions) ; autre - 
ment dit, la repartition des vitesses est invariante par translation suivantz, soit : V(M) = u(y)x 
Comme le fluide est visqueux, les parois des plaques fixes imposent une condition limite a la vitesse : la vitesse 
du fluide qui touche une paroi doit etre egale a celle de la paroi : c’est la condition d'adherence. 

Lorsque le fluide est visqueux, il faut rajouter dans le 2 e membre de l’equation du mouvement (c’est l’eq. 
d’Euler si le fluide est parfait) les forces de frottement d’origine visqueuse. Si le fluide est newtonien, sa visco- 
site est constante pour une temperature donnee. Les efforts visqueux sont alors proportionnels aux gradients de 
vitesses. La formulation de ces efforts visqueux est simple dans le cas de 1’ ecoulement entre deux plaques car le 
champ de vitesse est unidirectionnel, V(M) =u(M)x . 




Cas ou les deux plaques sont fixes. 

On a comme conditions limites : 



y = 0 => u = 0 



y = a => u = 0 



(XIII-2) 



on cherche une vitesse V dont les composantes sont telles que : 



u = u(y) v = w= 0 (XIII-3) 

Apres calculs, on trouve pour la composante u : 



u{y) = 



a 2 dp g 

2p dx 




y_ 

a 



(XIII-4) 



dp / dx est appele gradient de pression motrice de 1’ ecoulement ; c’est une chute de pression motrice propor- 

tionnelle a la longueur parcourue : la chute de pression divisee par la longueur parcourue est constante car 
1’ecoulement est etabli (la pression motrice decroit lineairement, voir la figure ci-dessous). Ainsi, entre deux 
sections distantes d’une longueur L on a : 



dP g _P g 2~P g 2_ A P g 



(XIII-5) 



dx x 2 —x l L 

Pour determiner la derivee de la pression motrice par rapport a x il n’est done pas obligatoire de se ramener au 
voisinage d’un point. La distinction des deux cas ou dp dx = Cte et dp g / dx = f (x) sont illustres ci-dessous. 
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Entre les deux plaques, on a une distribution de contraintes visqueuses donnee par : 



1 dp 



12p dx 



g 3 

a 



(XIII-6) 



Le debit volumique est L integrate de u sur une section droite ; 



q v = \ 0 U ^ dy \ dz = ~ 1^7 -fa- 1 dz (XIII-7) 

z z 

On choisit d’exprimer ici le debit par unite de largeur (la largeur est perpendiculaire a la figure, soit z) ; ce debit 
unitaire est donne par : 



c Jv 

l dz 



dv 




a 3 dp g _ a 3 &Pg 
12p dx 12p L 



(XIII-8) 



La vitesse moyenne de l’ecoulement est obtenue a partir du debit et de l’aire A d’une section de passage de refe- 
rence S perpendiculaire a l’ecoulement: 




y _ dv _ a g 

q ~ A ~ 12(i L 



(XIII-9) 



Le point ou la vitesse est maximale est situe a : 

y = j => ^=\v q (xm- io) 



La vitesse maximale est ici 1,5 fois plus grande que la vitesse moyenne. 



XIII-4 Ecoulement laminaire developpe etabli dans une conduite de section circulaire 



La « longueur d’entree », L E , est la longueur de tuyau entre l’entree du tube et la position ou fecoulement lami- 
naire est completement developpe et etabli. Cette longueur est donnee par une formule empirique : 



L e = 0,058 ReD 



(XIII- 11) 
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ou D = 2R est le diametre interieur du tube et re, le nombre de REYNOLDS de L ecoulement, defini a partir de 
la vitesse moyenne V q du diametre D, de la masse volumique p et de la viscosite dynamique p du fluide : 



Re = 



P V q D 
P 



(XIII-12) 



La vitesse de P ecoulement etabli est ici aussi de la forme : 

V = ux (XIII- 13) 

oil Ton peut montrer que u ne depend que du rayon r dans une section droite ( 0 < r < R ). 
La vitesse verifie la condition d’adherence sur la paroi : 

r = R => u = 0 (XIII- 14) 



Le profil des vitesses dans la section droite est : 

1 dp g 



ii (r) = 



4p dx 



V-K 1 ) 



(XIII- 15) 



Le gradient de pression motrice est ici aussi une constante pour un ecoulement donne ; on a done, entre deux 
sections distantes de L : 



Oq 

1 


A P S 


dx 


L 



(XIII-16) 



Le profil de vitesse peut done s’ecrire : 



u(r) = 



P* 

4p L 



R- 



( r 2 A 

1 R- 



(XIII-17) 



la distribution de contrainte est donnee par : 

du r OP x 



™ ^ dr 2 dx 



(XIII- 18) 

Le debit en volume verifie la formule celebre ci-dessous, connue sous le nom de « loi de POISEUILLE » : 



dv ~ ' 



n R a Ap 



8p L 



(XIII-19) 



On remarque que la loi de Poiseuille donne une relation de proportionnalite entre : le debit et le gradient de pres- 
sion motrice. Si on appelle «conductance» le terme : 




(XIII-20) 



(XIII-21) 



le gradient de pression motrice de P ecoulement, on a la relation : 
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(XIII-22) 







Qv 


q v = ao 


ou : 


a = — 






a 



On appelle conductance le terme a . 

La vitesse moyenne de L ecoulement est : 




(XIII-23) 



Le point ou la vitesse est maximale se situe sur l’axe du tube, soit en r = 0 . lorsque du / dr = 0 : 



r = 0 




dr 



u = u 



max 



R 2 A p g 

4q L 



= 2F 



(XIII-24) 



Dans un ecoulement laminaire etabli dans une conduite de section circulaire, la vitesse maximale est deux fois 
plus grande que la vitesse moyenne. 



XIV- Ecoulement dans les conduites en charge. 

XIV-1 Profils des vitesses dans les tuyaux 



Dans un ecoulement laminaire completement developpe , le profil des vitesses verifie la loi de Poiseuille : 

(XIV-1) 



1 A Ps 



“ (r)= tir* 2 






R- 



ou : 



u r " 



V„ 



R 1 



la vitesse est maximale au centre et nulle sur la paroi. 



Dans un ecoulement turbulent stationnaire, le profil des vitesses peut etre represente par la formule empirique : 




(XIV-2) 



oil l’exposant n est fonction du nombre de REYNOLDS que Ton a deja defini plus haut. 

Cette formule empirique est appelee « loi puissance ». On peut preciser quelques valeurs de l’exposant n : 
Re = 4xl0 3 , n = 6 

Re = 1,1 x 10 s , n = l (XIV-3) 

Re = 3,2 x 10 6 , « = 10 



Le rapport entre la vitesse moyenne et la vitesse au centre du tube est donne par : 



r\ 2 

u _ 2n 
T q ~ (77 + 1X277 + 1) 



(XIV-4) 



On voit que si n augmente, ce rapport augmente. 

Dans la figure qui suit, on a represente differents profils de vitesses pour fecoulement laminaire et un ecoule- 
ment turbulent dont le nombre de Reynolds est de 4000. 
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Conduite 



r/R ‘ 




u \ 

max 




Laminaire (Poiseuille) 






1 


Turbulent (loi puissance) 



XIV-2 Repartition des contraintes de cisaillement (ou frottement) dans les ecoulements 
etablis (cas des conduites). 

On fait les hypotheses suivantes : 

-Ecoulement horizontal. La pente du tube est nulle, done : dZ / dx - 0 . La pesanteur n’intervient pas dans le 
mouvement du liquide : dp g / dx = d(p + p gZ) / dx = dp/ dx + p gdZ / dx = dp / dx . 

-Ecoulement stationnaire 
-Fluide incompressible 

-Ecoulement completement developpe (etabli) 

L’ecoulement etant stationnaire la projection des forces exterieures au volume de controle suivant l’axe x est 
nulle : 



Conduite 




y f - x = o 



^ F ■ x = {p x - p 2 )nr 2 + x r 2nrdx = 



, , , , dp dx 



nr 2 - 



dp dx 
^*52 



nr 2 + xlnr dx = 0 



(XIV-5) 



D’ou : 






r dp 
2 dx 



(XIV-6) 



En particulier, sur la paroi du tube ( r = R) , 



_ R dp 
rx ~~^2~dx 



(XIV- 7) 



Mecanique_des_fluides_2A_MP.doc 



27 



21/11/2007 






XIV-3 Bilan d’energie dans un tube : PERTE DE PRESSION TOTALE et PERTE DE 
CHARGE 



Considerons une conduite horizontale comportant un changement de direction ou « coude ». 

Si on applique le premier principe de la thermodynamique au volume de controle D limite par la surface S, on a 



P -l- P P _1_ P — JT 

x S ' ± frottement ' x autre 



= ^|p(M)e(M)c?a)+|je(M*) + j p(M*)v(M*)] P V(M*)-cIa\ (XIV-8) 
oil le premier membre represente les puissances « produites » (au sens algebrique), et : 

(XIV-9) 



e = u + ^V 2 + gZ 



u est 1’ energie mecanique par unite de masse 



Si on fait les hypotheses suivantes : 



p - p - p =o 

S frottement autre 



(i) 

(ii) ecoulement stationnaire de fluide incompressible 

(iii) pression et energie interne uniforme dans le volume de controle, 

le bilan d’energie precedent se reduit a : 



1 , 
p l +pgZ\+-pa l v ql 



1 2 

p 2 + pgZ 2 + 2 P a 2 V ,2 



t \ 8 Q 

= P(«1 -w 2 j-p 



8 m 



(XIV- 10) 



ou : 



a est le coefficient d’energie cinetique 
V est la vitesse moyenne 



1 2 
fl+P?2|+yP a /,i 

1 9 

p 2 + pgZ 2 +-pa 2 l / „ 2 

5 Q 



est l’energie mecanique totale par unite de volume dans la section SI 
est 1’ energie mecanique totale par unite de volume dans la section S2 



p(w, — m 2 )— Pg^-est la difference d’energie mecanique totale entre les sections SI et S2. 
On appelle perte de pression totale le terme : 



sip) 

AP?i_ 2 = p(w| -u 2 )-p 



8 m 



(XIV- 11) 



et perte de charge le terme equivalent exprime en hauteur de fluide : 



A H x 2 = APt Al_ 

Pg 



(XIV- 12) 



On peut done ecrire : 



1 2 
p l +pgZ l +-pa l V qi 



1 2 
p 2 +pgZ 2 +jpa 2 V q2 



= A Pt x _ 



(XIV-13) 
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On remarquera que l’energie mecanique totale dans une section est ecrite a l’aide de la vitesse moyenne V q dans 

cette section alors meme que le profil des vitesses n’est pas « plat » (voir les fonnules donnant les vitesses). Ceci 
n’est possible que grace a 1’ introduction du coefficient d’energie cinetique a qui « corrige » la valeur moyenne 
de l’energie cinetique. 11 n’y a pas lieu, par contre, d’introduire de coefficient de pression pour le terme p + pgZ 
puisque celui est constant dans le plan de la section droite. 

XIV-4 Perte de charge reguliere et perte de charge singuliere 

La perte d’energie mecanique totale dans un tuyau se decompose en deux parties : la perte de charge reguliere 
qui se produit tout le long de fecoulement et les pertes de charge singulieres qui sont localisees en des endroits 
precis du conduit (coude, valve, raccords, etc..) 

Pour tenir compte de ceci, on ecrit : 



k-Ph-2 - (a Pf|_ 2 ) r4g + (a Pt x 2 ) sing 



(XIV-14) 



XV - Calcul des pertes de charge. 

XV-1 Coefficient de perte de charge reguliere 

Quel que soit le regime de fecoulement (laminaire ou turbulent developpe), la perte de charge peut etre mise 
sous la forme suivante : 






(XIV-15) 



C’est la formule dite de DARCY- WEISBACH. 

Cette formule est obtenue par analyse dimensionnelle, en appliquant le theoreme de VASHY-BUCKINGHAM. 
Elle donne la chute de pression totale se produisant entre deux sections distantes de la longueur L. et en ecoule- 
ment completement developpe. 

D h est le diametre hydraulique ; il est defini a partir de l’aire A de la section droite et du perimetre mouille 
X Par: 




(XIV-16) 



On appelle coefficient de perte de charge reguliere le coefficient sans dimension A ; c’est une fonction du 

nombre de Reynolds et de la rugosite e interieure du tube (si l’ecoulement est turbulent): 



A = f(Re,e) = f x Re k , 




(XIV-17) 



Cette fonction de Re est notamment donnee dans le diagramme de MOODY. On peut citer aussi les travaux de 
NIKURADSE concernant la determination de cette fonction. 

La determination de cette fonction consiste a trouver le coefficient / et l’exposant k du nombre de Reynolds. 

On appelle coefficient de forme le coefficient sans dimension f 11 ne depend, pour un regime d’ecoulement 
donne, que de la forme de la section droite du tube. 
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Exercice : Etablir la formule de Darcy-Weibach pour un ecoulement etabli en supposant que la chute de pression totale A P t 
peut-etre ecrite comme une fonction de : D h , p . L , V , pour un tube a paroi lisse, fonction telle que : 

AP l =f(D h , l i, P ,L,V q ) = CD a y p c L“ (1) 

ou C est une constante sans dimension. 



II s’agit done de trouver les exposants a,b, c, d, e, de telle sorte que la dimension de C D“ t \l h p c L d F e soit identique a 
celle de A P t , soit : [a.P,] = ML *T 2 . 

Les dimensions de chacun des tennes influenqant la perte de charge sont : 

[Aj=L; [p] = ML- 1 T- 1 ; [p] = ML- 3 ; [fJ=LT 1 ; [,i] = ML^T" 1 . 

En exprimant que la formule ( 1 ) est homogene du point de vue dimensionnel : 

ML 'T 2 = (l)“ (ML 'T 2 ) 4 (ML 3 ) c (L) (LT- 1 ) e . 

Experimentalement on montre que dans un ecoulement etabli, la perte de charge est proportionnelle a la longueur : on peut 
done ecrire : d = 1 . 

En developpant puis en identifiant respectivement les exposants de chaque grandeur (Masse, Longueur, Temps), on a : 
b + c = 1 a = e — 3 

a-b-3c-e = -2 => <b = —e — 2 

—b—e=— 2 c=e— 1 

a, b, c, sont parametres ici par e. D’ou : 

F e . 

<1 

Faisons apparaitre dans la perte de charge le nombre dc REYNOLDS de l’ecoulement ainsi que l’energie cinetique 
moyenne p F 2 / 2 : 



a p t = c rr 3 p - e+2 p 2 l 



A P t =C D"- 3 p- e+2 p e - 2 L 



v e 



= 2 C 



e-2ye- 2 j^e-2 

q h 



LD^-p V ;=2CR^j-\pV; 



En posant : 



f = 2C 
k = e — 2 
A = /'Re* 

on obtient : 



AP, = A — — pF 

D „ 2 



Coefficient de forme 

Exposant dependant du regime d’ecoulement 
Coefficient de perte de charge 



XV-1-1 Cas de l’ecoulement laminaire dans une conduite de section circulaire cons- 
tante : loi de Poiseuille 



L’ecoulement de Poiseuille possede une solution exacte ; il est done possible d’obtenir theoriquement le coeffi- 
cient de forme directement a partir du proftl des vitesses u. 

En effet, 





(XV- 1) 



L’integrale de u figure dans la vitesse moyenne, puisque : 
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V q = ~7 f u(M)dA 
* A J MsS 



(XV-2) 



Les courbes experimentales montrent que si l’ecoulement est laminaire developpe, 
k = -1 (XV-3) 

On a done, en se rappelant que V q = q v / A , avec A = n R 2 et x = 27 iR : 



f 





pV q D h APt m 2 



AP/ 32A 3 
L PX 2 9v 



APt 32n 2 R 6 
1 nr <lv 



On exprime le debit grace la loi de Poiseuille en utilisant : 



q v = ao 



APt nR 4 
L 8^ 



D’ou : 



/ = 



APt 32nR 6 



L 



\\An 2 R 2 



APt n R 4 
L 8(j. 



Apres simplifications, on obtient la valeur entiere : 



(XV-5) 



(XV-6) 



(XV-4) 



/ = 64 



(XV-7) 



On notera bien que e’est une valeur exacte (et non deduites de mesures). 

L’ecoulement entre deux plaques fixes paralleles possede aussi une solution exacte (ecoulement de Couette) ; 
son coefficient de forme vaut : / = 96 . 



XV-1-2 Cas des conduites cylindriques de section non circulaires 

En ecoulement laminaire, et lorsque la section n’est pas circulaire, (cas des echangeurs de chaleur...), on peut 
calculer le coefficient de forme a condition que Ton ait au prealable determine le profil des vitesses u dans la 
section droite. 



Si on deflni les variables sans dimension suivantes : ( p viscosite dynamique du fluide, a est le gradient de 
pression motrice, a est une longueur caracteristique de la section droite, par exemple la longueur du cote si la 
section est rectangulaire) 

u* = " 2 u (vitesse normalisee) 

a a 2 

x* = x/ a , y*= y / a , z* = z / a (coordonnees normalisees) 

S devient S * (section en coordonnees normalisees) 



l’equation du mouvement permanent du fluide suivant l’axe du tube se ramene a une equation de Poisson : 



Au * = -1 



dans la section de passage, 



: = 0 , sur le bord de la section de passage S* (condition d’adherence sur la paroi) 
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La resolution de cette equation de Poisson est possible par des methodes numeriques (differences finies ou ele- 
ments finis si le contour est «complexe»). Le profil des vitesses u* etant connu, on peut trouver par integra- 
tion la conductance (en valeur adimensionnelle) : 



cr* — — — g = \u*{M)dy*dz* 

n J 



((XV-8) 

D’ou la valeur du coefficient de forme (en notant A * l’aire de S* et x* = X / a le perimetre normalise de S*) 
/ - ' 7 ' (XV-9) 



c*x * 2 



XV-1-3 Cas de l’ecoulement turbulent etabli. 

Pour les ecoulements turbulents on doit utiliser des resultats experimentaux pour determiner le coefficient de 
perte de charge A en function du nombre de REYNOLDS Re. Ces resultats sont presentes dans un diagramme 
dit de MOODY, tel que celui de la figure ci-dessous qui s’applique aux conduites circulaires*. 




La « rugosite » est notee s, c’est le rapport entre le diametre interieur et la taille moyenne, t, des asperites : 
e = D / 1 . C’est un parametre important pour Lecoulement turbulent. 



Les differentes courbes correspondent aux formules suivantes : 



1 ^ Loi de POISEUILLE : 

2 — > Loi de BLASIUS, conduites lisses : 

3 ^ Formule de KARMAN-PRANDTL : 

4 ^ Formule de KARMAN-PRANDTL : 



. _ 64 . 
Re ’ 



A = 



0,316 

Re 1 ' 4 



-^= = 2 log 10 [Re 



Va] - 0,8 






= 2 log! 



D 



+ 1,14 



* Les resultats traces ici sont ceux de NIKURADSE. 
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En resume, si on veut calculer les pertes de charge regulieres d’un ecoulement en regime turbulent, on doit sui- 
vre la demarche suivante : 

(a) Evaluer le nombre de Reynolds par la formule : Re = p V D h / p ; en regime turbulent. Re » 2000 . 

(b) Evaluer la rugosite interieure s de la conduite 

(c) En deduire le coefficient de perte de charge A grace au diagramme de MOODY ou a des formules empiri- 
ques fournies dans les ouvrages d’hydraulique. 

(d) Calculer le perte de charge, soit en unite SI de pression (Pa), soit en hauteur equivalent de liquide, en appli- 

L 1 , 

quant : A Pt = A— — pV q 

Si la conduite est lisse, (rugosite nulle) on se rappellera que la loi de BLASIUS est applicable ( / = 0,316 et 
k = -1/4) pour le regime turbulent, soit : 



A = 0,316 x Re~ l/4 



(XV-10) 



XV-2 Coefficient de perte de charge singuliere 



La perte de charge singuliere est exprimee en fraction de l’energie cinetique du fluide a P entree de la singularite. 
On utilise la formule suivante si cette perte singuliere est exprimee en unite de pression : 



4 



Homogene a une pression ( ML *T 2 , Pa en unite SI) (XV-11) 

ou la formule suivante si cette perte singuliere est exprimee en hauteur equivalent de fluide : 

Homogene a une longueur (L, metre en unite SI) (XV-12) 






Le coefficient K est appele coefficient de perte charge singuliere ; il est sans dimension. Ce coefficient est par- 
fois note t, . 

Le probleme pour Pingenieur est alors de savoir determiner le coefficient K. 

De nombreuses etudes experimentales ont ete realisees et leurs resultats sont publies dans des bases de donnees 
pour Phydraulique. On peut ainsi obtenir pour la plupart des singularites usuelles (coudes, reductions, vannes, 
raccords, etc...) la valeurs du coefficient K en fonction des caracteristiques geometriques des singularites et du 
nombre de Reynolds Re . 

11 est toujours possible de determiner soi-meme K. Pour cela il faut mesurer simultanement la perte de charge 
entre l’entree et la sortie de la singularite et la vitesse moyenne en entree (ou le debit q v ), puis tracer le graphe 

de A Pt sing en fonction de p / 2 ; le coefficient cherche est alors la pente du nuage de points. 





j pV q 2 (Pa) 
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11 faut noter que pour le meme nombre de Reynolds, des singularity « homothetiques » auront le meme coeffi- 
cient K. 

Par exemple, pour le meme nombre de Reynolds Re, deux coudes « droits » (angle de 90°) ayant le meme rap- 
port rayon de courbure/diametre (soit ) auront le meme coefficient K. 




XVI - Quelques elements sur la theorie de la couche limite 

XVI-1 Introduction 

Les ecoulements de couche limite sont caracteristiques des fluides visqueux en mouvement relatifs pres d’une 
paroi solide. 

La couche limite est la region du fluide en mouvement proche d’une paroi solide. Cette zone est particulierement 
sensible aux effets de la viscosite. 

Une zone laminaire debute au bord d’attaque (x = 0) et son epaisseur augmente a mesure que l’ecoulement 
avance le long de la paroi. 

11 se cree alors une zone de transition laminaire-turbulent apres laquelle / ’epaisseur de la couche limite aug- 
mente. 

11 existe dans la zone turbulente une sous couche visqueuse dans laquelle les effets visqueux (lies aux 
frottements entre couches de fluides animees de vitesses differentes) sont predominants devant les effets 
d’inertie (lies a la mise en mouvement des masses fluides). 

On represente le developpement de la couche limite sur une plaque plane dans la figure ci-dessous : 



8(x) epaisseur de la couche limite 




Zone laminaire Zone turbulente 



On appelle epaisseur de la couche limite la distance 5 entre le fluide et paroi a partir de laquelle la vitesse du 
fluide est egale a 99% la vitesse U 0 du fluide libe (non perturbe par la paroi). 
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On appelle epaisseur de deplacement 8 * la longueur du deplacement qu’il faudrait appliquer a la paroi solide 
pour conserver le meme debit en masse a un ecoulement de fluide parfait hypothetique. 

On a, en fluide incompressible : 



8 * = 



-f 



1 - 



u(x M ,y) 

U n 



dy 



(XVI- 1) 



XVI-2 Equation integrate 

Admettant que la pression motrice s’identifie a la pression statique p dans la section droite et que de plus p est 
constante dans cette section, la resultante des forces suivant la direction x s’ecrit : 

dF x = a(x) - a(x + dx) + ct( 8) - x p dx (XVI-2) 

Frontier? de la couche limite Volume de controle 

n 

o 
s= 
o 
cr 

CD 

| 

r-t- 

CD 



Plaque 
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Ou x est la contrainte de frottement sur la paroi. 

En considerant une unite de longueur perpendiculairement a la figure, on a comme force : 



a(x) = px8, <j(x + dx) = 



p{x) + dx 
dx 



(5 + dd) 



o(8) = 



p(x)+ TTx dx 



x d8 



dF x =px8 



P(x) + \t dx 



: (8 + 1/8)- 



P(X) + lt dX 



d 8 - x p dx 



(XVI-3) 



Si Ton combine le flux de quantite de mouvement suivant x avec l’equation de conservation de la masse pour le 
volume de controle choisi, on obtient la forme generate de 1’ equation integrate de quantite de mouvement : 



~ dp , d 

~^~df~ X P dX ~~df 


f 5 2 
P U dy 

Jo 


-U — 
0 dx 


| p u dy 



(XVI-4) 



Si Ton veut evaluer l’epaisseur de couche limite pour un ecoulement laminaire de fluide incompressible au des- 
sus d’une plaque plane, il faut : 

(i) Determiner le gradient de pression dp / dx en dehors de la couche limite (par une etude en ecoulement 
potentiel) ; 

(ii) Se donner une expression analytique plausible pour le profil de vitesses u(y) dans la couche limite ; 

(iii) Utiliser comme formule pour la contrainte visqueuse (frottement), la loi de Newton : 



du 




XVI-3 Ecoulement dans la couche limite d’une plaque plane 

On prend les hypotheses suivantes : 



L’ ecoulement exterieur possede une vitesse constante : 
La pression est constante : 

Le fluide est incompressible : 



U 0 = Cte ; 




p = Cte . 
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L’equation integrate de quantite de mouvement se « simplifie » et conduit a : 




Soit la variable sans dimension definie en x donne : 

y 

r| = — => dy = 8dx\ 

o 

On obtient : 




(XVI-5) 



C’est l’equation integrate des quantities de mouvement pour un eeoulement sans gradient de pression ; 

cette expression de la contrainte de frottement est restreinte au cas suivant, rappelons-le : 

(i) Eeoulement stationnaire sans forces de pesanteur (le fluide est dit « non-pesant ») ; 

(ii) Eeoulement incompressible bidimensionnel ; 

(iii) Eeoulement le long d’une plaque plane. 



XVI-4 Equation integrate pour un eeoulement sans gradient de pression ( dp / dx = 0 ) 

A 1 Eeoulement laminaire le long d'une plaque plane 

On se donne la forme du profil de vitesses : 

u(y) = A + By + Cy 2 

c’est un polynome de degre 2 en y. 11 doit verifier les conditions limites : adherence sur la plaque (l ere condition), 
raccordement avec l’ecoulement exterieur de vitesse U 0 (2 emc et 3 eme condition) : 



II 

o 


=> 


u = 0 


y = § 


=> 


u = U 0 


y = 8 


=> 


du 

dy 



(XVI-6) 



En injectant l’expression polynomiale ci-dessus dans la formule de x p vue plus haut et en appliquant la loi de 
Newton pour le membre de gauche, on obtient : 




(XVI-7) 



Par definition, le coefficient de frottement est donne par le rapport entre la contrainte de frottement et l’energie 
cinetique (ou pression dynamique) de 1’ eeoulement exterieur : 



C/ = 




(XVI-8) 



Dans le cas present, on peut montrer que : 
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(XVI-9) 



0,730 
1 ^Re(x) 



Ou Re(x ) est le nombre de REYNOLDS a l’abscisse x : 

Re = (XVI- 10) 

P 



L’epaisseur de la couche limite est donnee par : 



§ _ 5,48 
* jRe(x) 



L’epaisseur de deplacement est donnee par : 



5* 1,835 

x 4 R e(x) 



(XVI-11) 



(XVI- 12) 



B / Ecoulement turbulent le long d'une plaque plane 

On se donne la forme du profil de vitesse : 



u 

U~ n 



y 

v§y 



(XVI- 13) 



la contrainte de frottement, determinee experimentalement (PRANDTL), verifie : 



x= 0,0225 p Ul 



pU 0 8 



(XVI- 14) 



L’epaisseur de la couche limite est donnee par : 



S _ 0,37 

* ~ Re(x) 1 ' 5 



(XVI- 15) 



Le coefficient de frottement sur la paroi est donne par : 



C, = 



0,0577 



Re(x) 



— , lorsque : 5x10 <Re(x)<10 



(XVI- 16) 



XVI-5 Coefficients de frottement des plaques lisses en ecoulement turbulent 

Pour des nombres de Reynolds superieurs a 10 5 , le coefficient de frottement C , est exprime par une formule 
deduite de resultats experimentaux : 

C f = 0, 074 Re L ~ h ' \ lorsque : 5x10 s <Re L <10 7 (XVI- 17) 

C f = °’ 455 . 2 58 , lorsque : Re L > 10 7 (d’apres H. SCHLICHTING) (XVI- 18) 

(log Re L ) 
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ou 



Re L = 



?U 0 L 

d 



est le nombre de Reynolds defini a partir de la longueur de la plaque. 



La figure ci-dessous represente le coefficient de frottement C f du fluide leger sur la paroi d’une plaque plane 
lisse en fonction du nombre de Reynolds Re L . 




1 E5 2 5 1 E6 1 E7 1 E8 1 E9 



L 

Courbe 1 : Ecoulement laminaire. 

Courbe 2 : Ecoulement de transition laminaire-turbulent 
Courbe 3 : Ecoulement turbulent. 
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